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ABSTRACT 

Under the unified model for active galactic nuclei (AGNs), narrow-line (Type 2) AGNs are, in fact, 
broad-line (Type 1) AGNs but each with a heavily obscured accretion disk. We would therefore ex- 
pect the optical continuum emission from Type 2 AGN to be composed mainly of stellar light and 
non-variable on the time-scales of months to years. In this work we probe the spectroscopic vari- 
ability of galaxies and narrow-line AGNs using the multi-epoch data in the Sloan Digital Sky Survey 
(SDSS) Data Release 6. The sample contains 18,435 sources for which there exist pairs of spectro- 
scopic observations (with a maximum separation in time of ~ 700 days) covering a wavelength range of 
3900 — 8900 A. To obtain a reliable repeatability measurement between each spectral pair, we consider 
a number of techniques for spectrophotometric calibration resulting in an improved spectrophotomet- 
ric calibration of a factor of two. From these data we find no obvious continuum and emission-line 
variability in the narrow-line AGNs on average - the spectroscopic variability of the continuum is 
0.07± 0.26 mag in the g band and, for the emission-line ratios log 10 ([N II] /Ha) and log 10 ([O III]/H/3), 
the variability is 0.02 ± 0.03 dex and 0.06 ± 0.08 dex, respectively. From the continuum variability 
measurement we set an upper limit on the ratio between the flux of varying spectral component, 
presumably related to AGN activities, and that of host galaxy to be ~ 30%. We provide the corre- 
sponding upper limits for other spectral classes, including those from the BPT diagram, eClass galaxy 
classification, stars and quasars. 

Subject headings: galaxies: general - techniques: spectroscopic 
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1. INTRODUCTION 

Active galactic nuclei (AGNs) are found to vary in the 
time-domain at many frequencies. The underlying physi- 
cal mechanism drivin g the variability seen in AGN spec - 
tra is unknown (e.g.. (Peterson] 120011: [Bono et al.ll2003ft . 
Several scenarios ha ve been proposed inclu ding accre- 
tion disk i nstability (iKawaguchi et al.lll998ft . supernova 
explosion (jAretxaga et alj 1 1997ft . microle nsing by stars 
in the intervening galaxy (|Hawkinsl 1 1 993ft . and varying 
ionizing source or/and varying optical depth of the ma- 
terial in the vicinity of light-generating regions (e.g., 
iTohline fe OsterbroclJT97alGoodrich|[T989ft . At present, 
models involving non-thermal emission from the jets or 
events occurring on the AGN accre tion disk are most 
favored (see, e.g.. IStalin et al.ll2004 for a summary of 
variability models in AGNs by radio loudness). 

Most variability studies, however, have been focused on 
Type 1 (broad-line) AGNs. While the emission lines of 
Type 2 (narrow-line) AGNs were found to vary in the X- 
rays on time-scales from hours to years (jGuainazzi et aTl 
ll99dlMueller et alJl2061lMateos et~aT1l2007ft . their opti- 
cal continua and emission lines are generally considered 
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to be non-variable. For the emissions from the narrow- 
line regions, the spatially low-density region (electron 
density n e ~ 1, 000 cm " 3 ) results in a l ong (~ 100 years) 
recombination time 7 (jPeterson et al.1 11995ft . exceeding 
the duration of a typical observation. For the opti- 
cal continuum, it is commonly considered to be com- 
posed mainly of the stellar light, in ac cord with the ex- 
pectation from the un ification model (jAntonuccil 119931 : 
lUrrv fc Padovani| [l995) in which Type 2 objects are ob- 
scured by a dusty torus along the line-of-sight. Domi- 
nated by the stellar light, the optical continuum in Type 
2 AGNs should therefore be temporally non-variable. 
There is, however, no strong observational evidence in 
the time-domain to either support or dispute these ideas. 

From a 4-year monitoring of the broad-band opti- 
cal variability of 35 Type 1 and 2 Seyfert galaxies, 
IWinkler et alj (|1992ft demonstrated that most galaxies 
in their sample were variable. Another promising work 
to probe Type 2 AGN variability in the optical is the 
intra- night variability of Seyfert 2 galaxies observed by 
iJand (]2001ft . where two o ut of the thre e objec ts varied 
by ~ 0.25 mag. Furth er. iTrippe et all (l2008ft found a 
Type 1.9 Seyfert galaxy (|Osterbrocklll 989) to change into 
a Type 2 Seyfert over a few years, and proposed the un- 
derlying cause to be a varying ionizing continuum. 

The UV-optical variability of broad-line AGNs is bet- 
ter established with an amplitude of ~ 10% (e.g., see 
discussions m IVanden Berk et all 12004ft . Wavelength 



7 The recombination time r e 
atom, so that r e = l/(n e a^) 



is referred to that of the hydrogen 
where n e is the electron density. 
The total recombination coefficient of hydrogen, a a, is equal to 
4.18 X 10~ 13 cm 3 s _1 at an electronic temperature T e = 10 4 K 
(Ostcrbrock 1989). 
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dependence of QSO variability has been studied by 
IWi hite et al.l <l2005fl. an d the C IV-emission dependence 
bv IWilhite et all (12006ft in t he Sloan Digital Sky Sur- 
vey ^SD^sTI^rkeiaIJl2003) ■ In these studies the au- 
thors have developed a method to remove wavelength- 
dependent systematics which may exist in the spectra, 
from a spectroscopic plate 8 at one epoch of observation 
to the same plate at another epoch. Correcting for this 
on a plate-to-plate basis they were able to select variable 
QSOs and to demonstrate a power-law-like QSO differ- 
ence spectrum; which in turn can be explained by the 
change in the accret ion rate in thermal accretion disk 
(jPerevra et al.l l2006f h QSO photometric variab i lity in 
the SDSS has been studied bv lVanden Berk et all (|2004f) 
for its dependency on physical paramete rs such as red- 
shift, luminosity and radio properties; by Ide Vries et al.l 
(2003) in which the authors used both the SDSS and 
the historical observations (up to ~ 50 years) and found 
the time-sca l e of t he va riability to be ~ 2 years; and by 
llvezic et al.l (|2004f) and iSesar et all (|2006ft in which the 
authors used the SDSS repeated imagings in combination 
with the Palomar Observatory Sky Survey (POSS) data 
to probe long term QSO variability (up to ~ 50 years), 
and constrained its time scale to be ~ 1 year in the rest- 
frame. 

Normal "inactive" galaxies (for our purpose defined to 
be galaxies that are not Type 1 or Type 2 AGN s, nor low- 
ioniz ation nuclear emission regions - (LINERS, [Heckman| 
1980)) arc not observed to be variable; and are generally 
not expected to be. For the continua of galaxies, the 
time scale of observation is negligible compared with that 
of stellar evolution (for solar-mass main sequence stars, 
lifetime w 10 10 years). For emission-line regions, that 
arise from low-density environments and H II regions, 
similar arguments can be made to those for narrow-line 
regions, with the understanding that the e lectron densit y 
is likely to be lower (n e ~ 1 - 1, 000 cm~ 3 . lSpitzei1ll978D . 

The goal of this work is, therefore, to probe spectral 
variability of galaxies and narrow-line AGNs using the 
SDSS multi-epoch observations. The SDSS data have the 
advantage of large sample size with the same spectropho- 
tometric reduction, which is critical to our work because 
we expect any variability present in galaxies to be small 
(compared with the amplitude of the QSO UV-optical 
variability, for example). Several techniques are used to 
refine the spe c trosco pic calibration, including those by 
Wilhi te~et al.l (|2005ft . The resultant galaxy variability 
measurement is compared with that of stars and QSOs, 
independently analyzed in this work. 

In <J2]we discuss the samples used. In ^3]we present the 
variability measure. In 21 we describe the refinement of 
the spectroscopic calibration. In SjS]we present the galaxy 
continuum variability, and the emission line variability of 
the variable candidates based on the continuum. In ^51 
we present results on variability analyses on stars and 
QSOs as sanity checks. In $7] we summarize the results 
and discuss future work. 

2. SAMPLES 

As part of the SDSS (|York et all l2000ft spectra are 
observed with fibers of 3" diameter (corresponding to 

8 A spectroscopic plate in the SDSS contains a set of spectra 
which are observed simultaneously. 



0.18 mm at the focal p lane for the 2.5 m, f/5 tele- 
scope, iGunn et "ail 120061 ). The spectral resolution R is 
~ 1800 in the observed frame 3800 - 9200 A. All sources 
are selected (or "targeted") from an initial i maging sur- 
vey us ing the SDSS camera as described in lGunn et al.l 
(11998ft with th e filter response curves as described in 
iFukugita etldl (11996ft . and usi ng th e imaging process- 
ing pipeline of iLupton et all (12001ft. The as tromet- 
ric calibration is described in iPier et aT] (|2003l ). The 
photometric system a nd calibrati o n are described in 
IFukugita etldl (1 1996ft. [Hogg et all J2001ft. ISmith et all 
(|2002ft and llvezic et al.l (|2004ft . lTucker et all ([2006ft . The 
targeting st rategy for the multi- object spectrograph is 
described in lBlanton et alj (|2003ft . 

We select our samples from t he Data Release 6 (DR6, 
lAdelman-McCarthv et aLll200l of the SDSS. The galaxy 
sam ple is constructed fr om the SDSS Main Galaxy sam- 
ple (jStrauss et al J 12002). by selecting spectra which are 
both spectroscopically classified as galaxies, and reject- 
ing spectra in which the redshift measurements failed 
or have not been made. The spectral classification (star, 
galaxy, QSO) is highly confident, with 98% of all the DR6 
spectra having consistent classification between the spec- 
trold pipeline (|SubbaRao et al.l l2002ft and the specBS 
pipeline (Schlegel et al. in prep.). From this sample, 
pairs of spectral observations are identified by match- 
ing sources in both RA and DEC to within 1", and in 
redshift to within 0.01. Spectra at both epochs are re- 
quired to have the same plate number but different date 
of observations, given in Modified Julian Date (MJD) 
by the SDSS. The selection results in 23,330 pairs. The 
subsequent sample selections (mainly for the purpose of 
increasing the signal-to-noise, S/N, of the spectra) will 
be described in 21 and $5] Table [1] lists the relevant spec- 
troscopic plates and the corresponding number of galaxy 
spectra. The sample covers ~ 3 — 700 days in the ob- 
served frame. 

The spectroscopically classified galaxies selected above 
do not contain any strong (equivalent width, EW >10 A, 
and > 3er detection in the height of the line) broad line 
full width at half ma ximum, FWHM >1000 km s" 1 ) 
SubbaRa o et alj [20021 ). The cutoff is below the tradi- 
tional line- width selection for broad lines (Type 1) galax- 
ies, w here the w i dth o f permitted line(s) > 2000 km s _1 
(e.g., iKomossal 12008ft . The sample therefore com- 
prise s galaxies of different eClass types (IConnollv et al.l 
1995), and narrow-line AGNs. A potential contami- 
nation to our gala xy sample would be the na rrow-line 
Seyfert 1 galaxies (jOsterbrock k, Poggel 119851 . and ref- 
erences therei n), expected to offer a direct view of the 
active nuclei (|Antonuccil Il993t lUrrv fc Padovanil 11995ft 
and can be variable in the optical (|Giannuzzo &; Stirpg 
119961 : iPeterson et al.|[2000l ). Each of these spectra shows 
a broad component in Balmer emission line with FWHM 
< 2000 km s _1 , and relatively weak [O III]A5008. The 
reasons behind this potential contamination are twofold 
- firstly, the spectrold pipeline fits to each line a single 
Gaussian, that may underestimate the line width of the 
broad component in the case of a composite line which 
is made up of both the broad and narrow components, 
in particular when the broad one is relatively weak. Sec- 
ondly, the FWHM of the broad component of both Ha 
and H/9 of a narrow-line Seyfert 1 can extend down to 
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~ 500 km s" 1 (Fig. 4 of lZhou et al]|2006j) . falling within 
the galaxy spectral classification in the SDSS. However, 
we expect the number contribution of these objects to 
the whole DR6 ga laxy+QSO s ample to be small, ~ 0.5%, 
with reference to IZhou et al.l (|2006l ) who used an upper 
limit of 2200 km s _1 in the broad component of Ha or H/3 
emissions to search for narrow-lin e Seyfert 1 in the SDSS. 
Similarly, the Seyfert 1.8 and 1.9 (|Osterbrockll989t ). both 
of which show strong narrow components (relative to the 
broad ones, if available) in H a and H/3, can be variable in 
the optical ([Goodrich! fl989). and may be present in our 
galaxy sample as well. We note that the current anal- 
ysis does not provide further separation of those from 
our emission-line galaxies, that calls for fitting of double 
Gaussian to each Balmer line, for example. The number 
contribution of these types to t he full galaxy samp le is 
expected to be small, - 0.02% flWang L Weil 120081) . In 
this work, the fraction of the classified Seyfert 2 + star- 
forming + composite galaxies to the full sample of our 
galaxy spectral pairs is ~ 14%; and that of all the eClass 
types, ~ 22% (calculated from Table gj). 

When constructing the stellar and QSO samples, the 
spectral pairs are selected similarly to the case of galax- 
ies, except that they are spectroscopically classified as 
star and QSO, respectively. The targeting selection for 
spectrosco pic observation of QS Os in the SDSS is de- 
scribed in Richard s et al.l (12002'). Each QSO spectrum 
classified (jSubbaRao et al.|[2002l ) has at least one strong 
line (see above for the galaxy classification) with FWHM 
>1000 km s _1 , or/and a Lyman alpha forest. Unlike 
the selection of the galaxies and the QSOs, no criterion 
is imposed on the quality of the redshift when select- 
ing the stars. The upper bou nd of their distances i s 
cz = 450 km s _1 , or z = 0.0015 (jSubbaRao et alJl2002T ). 
The samples contain 9,100 stellar pairs and 3,205 QSO 
pairs (0.08 < z < 3.51), respectively. 

The spectra are de-reddened against Galactic extinc- 
tion using the library written by Simon Krughoff 9 , w hich 
adopts the SFD dust m aps (fSchlegel et ail 119981 ) and 
the extinction curve by IQ'Donnelli (| 19941) . Following 
the SDSS convention the spectra are expressed in vac- 
uum wavelengths. Flux densities fx are expressed in 
10~ 17 ergs s _1 cm -2 A -1 . 

3. VARIABILITY AMPLITUDE: F var 

We adopt the dimensionless variability measure, 
F var , for iV-epoch r epeat e d observations (e.g., 
iRodriguez-Pascual et all 119971 ; iPetersonl |2001|) . It 
is a fractional root-mean-squared variability amplitude 
defined as 



f = — — m 

where a 2 is the variance of the flux, S 2 is the mean square 
uncertainty of the flux, and (/) is the mean flux 

1 N 

ff2 =]vE(/*-</» 2 > (2) 
»=i 

9 The library is available from the author upon request. 
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The numerator, proposed by iBonoli et al.l (| 19791) . 
was termed the squ are-root of the "excess variance" 
(jVaughan et al.l 120031) . in the sense that the flux uncer- 
tainty is subtracted in quadrature from the variability. 
In the calculation of the continuum F vaT the wavelength 
intervals encompassing emission-line center wavelengths 
± 280 km s _1 are excluded (using the vacuum wave- 
length values in the line list 10 adopted by the SDSS, and 
excluding the lines Ca K and H, US, the G band around 
4306 A, MgA5177, NaA5896, and the Ca II triplet that 
usually appear as absorptions). The mean square uncer- 
tainty for the spectrum at each epoch (5f for a given i) is 
calculated by quadrature summation of the pipeline flux 
uncertainty per pixel, over all valid pixels. It is therefore 
a variance in flux based on photon statistics. A valid 
pixel is defined to be a pixel that is not flagged as bad. 
The flags are listed in t|4.2l F yax is calculated in the ob- 
served frame with wavelengths 3900 — 8900 A for each 
object and for all spectral types. 

When calculating F val for narrow-line AGNs, since any 
possible host-galaxy contribution is not explicitly sub- 
tracted from an observed spectrum, the full observed 
spectrum is used in calculating the denominator (/) . As 
the numerator term "excess variance" does not contain 
the non-variable host-galaxy contribution, the resultant 
F var is therefore a lower limit on the AGN variability. 

If the excess variance a 2 — S 2 < 0, we assume the vari- 
ability is zero for the object, but continue to include it in 
the analysis. In this work, N = 2. We note that any non- 
zero -F V ar is referred to as a "variability" measurement, 
regardless of the origin being physically interesting or 
otherwise. 

4. REFINING SPECTROSCOPIC CALIBRATION 

The principal steps in the spectrophotometric calibra- 
tion of each spectrum in the SDSS are: a wavelength- 
dependent calibration using the observed standard stars 
on the same plate; and the tying of the absolute flux 
scale to the observed PSF magnitudes of st ars, on the 
same plate (jAdelman- McCarthy et all l2008l K The re- 
sultant average uncertainty in the spectrophotometry 
for DR6 is 7% at observed fra me wavelength 3800 A 
(jAdelman-M cCart hv et all 12008). As any variability in 
the sources from our sample is expected to be of low am- 
plitude, the calibration of the pipeline spectra must be 
improved using the approaches described below. 

4.1. Wavelength-dependent correction 

We adopt the method developed by IWilhite et al.1 
(2005) to remove systematics which may be present in 
each pair of galaxy spectra taken on separate MJDs. 
Flux-density corrections as a function of the observed- 
frame wavelength in the range of 3900 - 8900 A are cal- 
culated by a linear fit (tied to the origin) between the 
flux density (/a) of all the objects at epoch 2 (the later 
epoch) and those at epoch 1. This calculation results 
in a calibration spectrum for one plate. We adopt a 
slightly different methodology from lWilhite et ail (2005) 

10 IStoughton eTHl ||200t ). 
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in which stellar spectra in each plate were used. Our cali- 
bration spectrum per plate is calculated using all galaxies 
present in that plate, as we want to maximize the number 
of objects in generating the calibration spectrum, and 
avoid potential systematic effects when applying point 
source calibration to extended sources. The underlying 
assumption is that the majority of the galaxies observed 
in a given plate are non-variable in time across the two 
epochs. Even if the galaxies were in fact variable, the 
method is still applicable because the variable ampli- 
tude of several hundred galaxies is not expected to be 
in phase. An example calibration spectrum is plotted in 
Fig. [T] Typically the calibration spectra are smooth. To 
focus on low-order corrections to the galaxy spectra, and 
to remove no i se in the calibration spectrum, we follow 
Wilhi te~et al.l (|2005| ) and set the final calibration spec- 
trum per plate as its 5th-order polynomial least-squares 
fit. By visual inspection this functional form works well 
in removing higher-order spurious spectral features. 

When constructing the calibration spectrum we drop 
any plate in which there are less than 10 pairs of galaxy 
spectra. We then refine each and every spectrum present 
in a given plate by using the calibration spectrum, if 
available. 

4.2. Wavelength-independent correction 

After removing wavelength-dependent systematics, 
each spectral pair can, in principle, be subjected to a 
wavelength-independent normalization in their flux den- 
sities. In order to obtain an accurate variability measure- 
ment of each spectral pair, the total continuum fluxes of 
the two spectra at a chosen common wavelength region 
are set equal. Specifically, for each pair the galaxy spec- 
trum at the epoch of lower S/N is shifted in a wavelength- 
independent fashion according to 

/ A a°-S/N) (correc ted) - C band / A (lowS/N) (original) , 

(5) 

in which Cband is a constant for each spectral pair, cal- 
culated as 

n — W (highS/N) /it. (low S/N) (P \ 

(-/band — -TAGband /^Aeband , (.0) 

where -FAeband is the total flux in the restframe wave- 
length range (which is referred to as "rest-band" for con- 
venience). In this work the rest-band is chosen to be 
6450 A ± 520 km s _1 , as such no prominent emission 
lines are present in the galaxy spectra 11 . Only the galaxy 
spectra in the first epoch are re-calibrated because it is 
typically the existence of low S/N data that results in a 
second observation. The combined procedure of system- 
atics removal and restframe band scaling is referred to 
"SYS+BAND" in the following discussion. 

Given an initial sample (<j2l 23,330 objects), we im- 
pose S/N cuts using the rest-band measures (or the 
[O III]A5008 (vacuum) line, described later when we dis- 
cuss absolute flux calibration). In this case we reject any 
spectral pair if there are one or more pixels flagged in 
either or both spectra in the wavelength region of in- 
terest, and we require the flux to be larger than zero 

11 When refining QSO spectra which are located at a 
wide range of redshift, an observed frame wavelength range 
6450 A ± 520 km s" 1 is used. 



at both epochs. The flags under consideration are NO- 
PLUG, BADTRACE, BADFLAT, BADARC, MANY- 
BADCOL, MANYREJECT, LARGESHIFT, FULL- 
REJECT, SCATTERLIGHT, CROSSTALK, NOSKY, 
BRIGHTSKY, NODATA, COMBINEREJ, BADFLUX- 
FACTOR, BADSKYCHI and REDMONSTER. All of 
the spectral flux or S/N calculations in this paper are 
carried out with pixel-masking using these flags. After 
this S/N criteria, the number of objects in the galaxy 
sample is reduced to 18,435 (Table [1]). 

In Fig. [2] and [3] we show the restframe galaxy spectra 
in two epochs from the SDSS pipeline and after the cal- 
ibration refinement, respectively. Firstly, we note that 
the pipeline performs well, that the fractional difference 
between both epoch is about 10% at the shortest wave- 
length, which is sufficient for many spectral analyses in- 
volving emission or absorption lines. However, contin- 
uum residual in Fig. [2] is present, which induces a vari- 
ability (Fvar = 0.058). After the calibration refinement 
we are able to obtain a zero continuum variability. 

For comparing with the rest-band scaling previously 
discussed, w e also adopt the abso lute flux calibration 
described in iPeterson et al.l (|1998f ). In this procedure, 
the spectrum in the low S/N epoch in the duo is 
scaled wavelength-independently so that its total flux in 
[O III] A5008 (vacuum) is the same as that of the spec- 
trum in the high S/N epoch. The physical justification 
is that the [O III] A5008 line generally appears as a nar- 
row line arising from a low-densit y region, which is be - 
lieved to be non- variable in time (|Peterson et al.l fl995). 
We note that the average [O III] EW of galaxy spectra 
in the full DR6 catalog is a few A, so for some objects 
mainly the continuum flux is measured and is expected 
to be non-variable as well. The combined procedure of 
systematics removal and absolute [O III] calibration is 
called "SYS+OIII" in the following. The region of influ- 
ence is set to be 5008 A ± 280 km s _1 . 

When judging which method performs better in ana- 
lyzing spectral variability, we consider the better method 
to be the one which minimizes the mean, one standard 
deviation of the mean (o~ m ) and the one standard de- 
viation sample scatter (c) of the variability at a fixed 
spectral S/N. These statistics for the galaxy sample 
by using the SYS+BAND, SYS+OIII and other refine- 
ment methods are listed in Tabled] At S/N > 10 the 
SYS+BAND method produces smaller mean and scatter 
in F'var- The difference is not related to the difference 
in the width of the rest-band being used (520 km s _1 
vs 280 km s _1 ), as the SYS+BAND method is also re- 
performed using 6450 A ± 280 kins" 1 , the same width as 
in the SYS+OIII method. These statistics remain lower 
in value than those in the SYS+OIII even for this re- 
duced interval. The result suggests that the [O III]A5008 
line in some spectra may be impacted by the effect of 
seeing so that the total observed line flux is different for 
the two epochs. The other possible reasons for this dis- 
crepancy are: the total line flux is sampled at only a few 
wavelength bins in the medium resolution SDSS spectra, 
or the line may exhibit intrinsic variability. 

To further test the method SYS+BAND, we divide the 
galaxy sample into RED and BLUE galaxies by using 
the eClass classification (described in £)5.2.ip . The RED 
spectra are defined as objects with classification param- 



5 



eter <f>KL > 0°, and the BLUE galaxies, §kl < 0°. For 
each plate, the calibration spectrum is calculated using 
only the RED or BLUE galaxies, which in turn is used 
to re-calibrate galaxies of all spectral types on the same 
plate. The purpose behind this is to determine whether 
the refinement method SYS+BAND would be biased by 
the presence of potentially different galaxy spectral types 
in the plate. We find that using a specific spectral type 
does not cause the variability statistics to change (Ta- 
ble ©. We therefore adopt SYS+BAND as the calibra- 
tion method, and use all available galaxy spectral types 
when calculating the calibration spectrum throughout 
this work. The "before and after" of this approach is 
shown in Fig. 0] where we see that the scatter in F val 
vs. S/N of our galaxy sample is reduced after the refine- 
ment. Comparing the average variability of our galaxy 
sample using the SYS+BAND approach with that from 
the original SDSS spectral reductions (Table [2]) we find 
that the spectrophotometry is improved by a factor of 
two. 

4.3. Plate-to- Plate Difference 

There are plate-to-plate differences seen in the above 
definitions of variable and non-variable sources, that is, 
the best-fit coefficients in Eqn. [7] are in general plate 
dependent. In Fig.[5]we show (F var ) of the galaxy spectra 
vs. average spectral S/N, per plate. The average F val per 
plate is larger at lower overall spectral S/N, suggesting in 
those plates the variability measurement is likely affected 
by the presence of noise in the spectra. By dropping 
individual spectral pairs with spectral S/N < 10 at the 
second epoch, the average and the sample scatter in F mT 
are reduced to 0.010 and 0.003, respectively. 

5. GALAXY VARIABILITY 

5.1. Defining Variables 

Following IWilhite et ail £2005) we define the variable 
candidates by objects having 

i^ var (S/N) > b exp (y^j + b 2 . (7) 

The S/N of each spectrum is calculated in the sec- 
ond epoch and by using all of the valid pixels as de- 
fined in H4.2I The constants b ,bi,b 2 are determined by 
least-squares fitting the above exponential form to the 
binned data ((S/N), (F var ) + n cr(F var )), where cr(F var ) 
is one standard deviation of F var in each bin, and for 
n = 0, 1, 2, 3, 4. Typically 10 bins in S/N are used. 

The variable candidates (for galaxies, stars and QSOs) 
in this work are defined as objects with > 3cr detec- 
tion. We do not sub-divide the distribution for detec- 
tions greater than 4er, though in the future we plan to 
explore in more detail the tail of the variable distribu- 
tions. For comparison we have also defined variables as 
objects located in a fixed upper percentile of F var in the 
plot F var vs. S/N. The > 3a cut is equivalent to a 2 — 6% 
upper percentile for a given S/N (Table [3]). 

5.2. Continuum 

To explore any dependence of spectral variability or re- 
peatability on the spectral type of galaxies, we divide the 
galaxy sample using two different spectral classification 



schemes: eClass (|Connollv et al.lll995[ ) a nd the "BPT" 
diagrams or the Osterbrock diag rams (Baldwi n et al.l 
119811: IVeilleux fc Osterbrockl 119871 ). These two schemes 
are motivated by different spectral features. The eClass 
is a spectral type describing the steepness of the con- 
tinuum slope and other higher-order spectral features 
by using a linear combination of several eigenspectra 
for each galaxy spectrum. The eClass parameters, 4>kl 
and 8kl, are related to the coefficients in the eigenspec- 
tra expansion (called eigencoefHcients), which are cal- 
culated by the Spectrold pipeline in the SDSS. Specifi- 
cally, the angle <t>KL co rrelates with the Ha emission EW 
(jMadgwick et a l. 2003) which is an indicator of star for- 
mation rate. The angle 6k l discri minates galaxy spectra 
exhibiting post- starburst activities (|Connollv et al.lll995l : 
I Yip et al.ll200l . 

On the other hand, the BPT diagrams separate the 
narrow-line AGNs from galaxies such as the star-forming 
galaxies by emission line ratios which indicate the physi- 
cal conditions of gas inside the galaxies such as the elec- 
tron density and the mean level of ionization (| Osterbrockl 
1989). The objects in these two spectral classification 
schemes are not mutually exclusive. We hereafter present 
results on variability within both kinds of classification. 

5.2.1. eClass 

The spectral vari ability as a function o f S /N per mean 
eClass galaxy type (Connoll y et al.lfl995h is shown in Ta- 
ble |H and Table [5] for only the variable candidates. The 
mean eClass spectra are galaxy spectra where each is 
averaged over a non-trivial size of subspace of t he two 
eClass parameters 4>kl and 6kl (|Ym et al.H200l . Lim- 
ited by the number of spectral pairs, < 3 in types E or F, 
only the types A, B, C and D are considered in the vari- 
able candidates (Table E|) . The continuum of the mean 
galaxy spectrum is reddest in the type A and bluest in 
the type D. 

Firstly, we note that the average variability amplitude 
are consistent among different eClass types to within 1- 
2cr sample scatter. If we assume a priori that some pop- 
ulations of galaxies do not vary in time, then a spec- 
tral repeatability can be assigned to the population with 
the smallest mean and scatter in F waI , that is, the mean 
eClass type B. The spectral features of this class are sim- 
ilar to the elliptical gala xies in the atlas of nearby galax- 
ies bv iKennicuttl (jl992). The repeatability, taken to be 
the average F var , is 0.010 or 1.0% in the observed- frame 
wavelength range 3900 - 8900 A. 

5.2.2. Emission-line galaxies: narrow-line AGN vs. 
star-forming 

Emission-line galaxies are considered by using the BPT 
diagram, using the s tarburst theoretical modeling line by 
iKewlev et al.l (|2006|) . We consider the types: Seyfert 2 , 
star-forming, composite and LINERs (jHeckmanl 1 1 980h . 
We do not find any LINER in our galaxy spectra l pairs , 
using the classification criteria bv IKewlev et alj (2006, 
their Eqn. 15) which involved [O I]A6302. We hereby con- 
sider the diagram logi n (JO III] /R0) vs. logi ([N Ilj/Ha) 
(Fig. 5 of lBaldwin etalTll98l . 

In the first step of this analysis the EWs are taken 
from the SDSS reduction. Galaxies with any one of the 
four lines (H/3 at 4863 A, [O III]A5008, Ha at 6565 A, 
[N II]A6585) having restframe EW smaller than 1 A are 
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rejected (in the SDSS convention EW > for emission 
lines); and at least a 2a detection is required in each of 
the four lines. Restframe stellar absorption of the related 
hydrogen emission lines in the BPT diagram is corrected 
for by a constant i ncrement of 1.3 A (jHopkins et al.ll2~003l : 
iMiller et al.ll2003t l. Any possible aperture effect on the 
emission lines is neglected in this paper, because the 
primary interest of this work is the comparison of line 
strength between two epochs, which is done using the 
same aperture for each object in each epoch and will 
only be subject to the difference in seeing. To correct 
on an object-by-object basis would require an assump- 
tion of an average line strength as a function of radius 
for each galaxy. Although the true line-flux correction of 
individual objects should depend on the distribution of 
the gas and its projected angle on the sky relative to the 
observational aperture, under the typical sky conditions 
in the SDSS observatio n (the median PSF wi dth = 1.4" 
in the r band, Fig. 4 of [Abazaii an et aT1l2003f ) and given 
the diameter of the fiber, the seeing-induced aperture 
correction factor to [O III] is close to unity, based on a 
nearby (z = 0.009) Seyfert 1 with a well-resolved ( ~ 10") 
narrow-line region (Fig. 3 of lWanders etal] 11992^ . This 
means that seeing-induced differences in the [O III] flux 
between any two epochs should be small. Further, the 
SDSS 3" diameter fiber spectroscopy does not cause sub- 
stantial aperture bias in classifying g alaxy spectra based 
on emission lines ([Miller et al.N2003l) . 

In Table [4] we compare the variability amplitude be- 
tween the star-forming galaxies and the narrow-line 
AGNs in the full galaxy sample, and Table [5] for the 
variable candidates. It is a little surprising that the av- 
erage variability amplitude at S/N > 10 is 1.2% for the 
star-forming galaxies, slightly higher than that of the 
narrow-line AGNs, 0.9%. The difference is bigger than 
la m (i.e., the measurement of the F vaT difference is reli- 
able), but smaller than the la sample scatter, indicating 
the F val of both types are consistent with each other. 

In Table [6] we show the average absolute AB magni- 
tudes in g, r and i bands of variable candidates in the dim 
phase, and the sample-averaged magnitude difference be- 
tween the dim and bright phases, per spectral type. The 
magnitudes are calculated for each spectrum, by con- 
volving each restframe-shifted spectrum with the SDSS 
filter re sponse curves for ext ended sources at an airmass 
of 1. 3 dFukugita et aLlll996l) . and in the AB filter sys- 
tem (Ok e fc Gunnlll983D . The average difference in the 
g-band magnitudes ranges from 0.07— 0.10 mag depend- 
ing on the spectral type, but they are consistent among 
different spectral types to within la sample scatter of the 
magnitude difference, a(AM), a measure of the sum of 
the sample scatters in both the distributions of M(dim) 
and M(bright). For the narrow-line AGNs, no obvious 
continuum variability is found (e.g., 0.07 ± 0.26 mag in 
the g band). 

We note that the distribution of F var in the full sample 
is non-Gaussian. Using the > 3cr cut for the variable 
candidates, for a Gaussian distribution one would expect 
the ratio between the number of variable candidates and 
that of the whole sample to be 

1 - er/(3/V2) 

1 - erf(-oo) ' 1 ' 



or 0.13 %,where erf(z) is the error function, and 1 — 
er/(z) is the complementary error function of the form 
(2/0F) f 2 °° e-« dt. From Tabled we see that typically 
at various S/N the F var of the variable candidates span 
a longer tail (1.53 - 3.44 %). 

5.2.3. Structure Function of the Continuum 

The characteristic time scale of variability is usually 
studied by constructin g the structure f unctio n (SF, e.g., 
ISimonetti etaD (|1985h . iHughes et all (j!992f k and refer- 
ences therein), the variability amplitude F var as a func- 
tion of time lag. Although no obvious continuum vari- 
ability is found in our sample of galaxies on average, 
the SF is nonetheless constructed for the galaxy sam- 
ple. Through the SF analysis, an absence of a time scale 
would be a cross-check of the lack of variability, and the 
comparison between the SF of stars and of QSOs ($6[) 
with that of the galaxies can be carried out. 

The restframe SF, F var vs. Ar res t, of all galaxies in 
the galaxy sample is illustrated in Fig. [6j where Ar reS i = 
At/(1 + z), and z being the redshift of a given galaxy. 
The time lag At is the difference between the two MJDs 
of observation, MJD(epoch 2)-MJD(epoch 1). No ob- 
vious time-scale is seen in the above restframe SFs of 
galaxies for the — 500 days (restframe) sampled. A sim- 
ilar conclusion is obtained when the fractional change in 
flux (|A///|) is used instead of F var (Fig. [7]), where A/ 
is the difference in the flux between two epochs, and / is 
the average flux over the two epochs. 

An exponential function of the form 

F var = a (l-6e- AT "-/ r =) (9) 

is fitted to each structure function, where r s is the char- 
acteristic variability time scale. Typically a structure 
funct i on where b = 1 is considered (e.g., iBonoli et al.l 
119791 : iTrevese et al.lfl99l . Firstly we relax b to accom- 
modate the non-zero variability at zero time lag, obvious 
in Fig. [6l This can be attributed to the repeatability of 
the continuum variability measurement f £!5.2.1j) in eClass 
Type B, where the mean F var is 0.010 or 1.0%. The best- 
fit t s for the galaxies is ~ 10 years. Uncertainty of this 
time scale is large (~ 2 years, with a reduced \ 2 of 20) 
because the duration of the observation is much shorter 
than the best-fit r s . 

We then consider a typical b — 1 structure functional 
form at Eqn. [9[ By assuming F va-V at the zero time lag 
(-Fvar(O)) (the noise) is uncorrelated with the intrinsic 
-Fvar(the signal), one can obtain the intrinsic F vllI by sub- 
tracting F var (0) from the total F vax in a quadrature fash- 
ion. The resultant structure function is shown in Fig. [SJ 
and the best-fit variability time scales in Table [8j where 
-Fvar(O) is taken to be the minimum of the total F vav . 
The large best-fit \ 2 m stars and galaxies implies that 
the best-fit variability time scales are quite uncertain. A 
longer duration of observation would be ideal to improve 
the structure function. 

5.3. Emission Lines 

The BPT diagram is constructed for the variable can- 
didates to explore the narrow emission-line variabil- 
ity. Again the emission- line ratios log 10 ([N II] /Ha) vs. 
log 10 ([O III] /H/3) are considered. 

In measuring the EW of an emission line, we fit to a 
given continuum-subtracted spectrum a single Gaussian 
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function. The conti nuum is estimated by no n-negative 
lea st square fitt i ng (lLawson fc Hansonlll974h. provided 
bylDobos et all (|2006t) who used the iBruzual fc Charlotl 
(2003) stellar population 1 A-resolution models. In this 
way, the stellar absorption in the Balmer lines is implic- 
itly taken into account in the best-fit stellar continuum. 
In measuring the uncertainty of an EW, the following 
steps are performed. Firstly, the uncertainty in the lu- 
minosity densities of the stellar continuum of the spec- 
trum is calculated by using Sf^ = 8 fx, where S denotes 
the uncertainty of the continuum flux density f% and of 
the observed spectrum f\. The uncertainty in each EW 
is calculated to be ^ A S(f\ — f^)/f^d\, where A_r is 
the region of influence of the emission lines, taken to be 
around the line center ±280 km s -1 . Finally, the un- 
certainty in each of the line ratios log 10 ([N II] /Ha) and 
log 10 ([O III]/H/3) is propagated by using common error 
propagation formula. 

The line ratios between the dim and bright phases of 
sources are found to be located very closely on the BPT 
diagram, as shown in Fig. [9l The emission-line ratios 
and their uncertainties for various object types are given 
in Table [9l The uncertainty in the flux measurement is 
small, for example in the star-forming galaxies, which 
indicates a reliable line measurements. For both line ra- 
tios the difference between the dim and bright phases 
is typically comparable to the la sample scatter in the 
line-ratio difference (i.e., the last two columns of Ta- 
ble. [9|). This is true for all of the object types consid- 
ered, meaning we find no evidence of narrow emission- 
line-ratio variability which is above la level. When divid- 
ing the line-ratio measurements by the first and second 
epoch (Fig. [T0|) instead of the dim and bright phases, we 
also see no difference between the two that is larger than 
la. Further, the line ratios of the variable candidates 
and the non-variable candidates (Fig. [TTj) do not show 
substantial difference in the BPT diagram, true for all of 
the emission-line galaxy types. 

The above analyses are also repeated with the intrin- 
sically de-reddened spectra using the best-fit color ex- 
cess (E(B-V)) (Yip et al. 2008 in prep.), by the non- 
negative least square fitting of stellar populati on models 
on the SDSS DR6 galaxies (jDobos et al.ll2006f ). Expect- 
edly, each of the average line ratios remains unchanged 
to the la level (not shown), because the central wave- 
lengths of the two lines for a given line ratio are close 
and is thus insensitive to dust reddening. 

6. COMPARISON WITH STAR AND QSO 

To evaluate further the performance of the calibration 
refinement, we refine the calibration in stellar and QSO 
spectra in a spectroscopic plate using calibration spec- 
trum constructed from galaxies of the same plate. The 
procedure is described in S]4l 

6.1. Variability in SDSS Stellar and QSO spectra 

Table [10] lists the statistics of F var in our sample of 
the SDSS stellar spectra. The improvement in the re- 
peatability after the SYS+BAND refinement is a factor 
of three. The average F val of stars is found to be slightly 
higher than that of the galaxies. 

Table QT] lists the statistics of F var in our sample of 
the SDSS QSO spectra. The QSOs are located within a 
redshift range of 0.08 - 3.51 with an average of 1.16. 



Upon the SF analysis the characteristic time scale of 
the QSOs is found to be 0.4 ± 1.0 years (Table El), 
which agrees wi th the time scales commonly found in the 
literature (e.g.. ICristiani et al.l ll99l llvezic et all 12004 
IVanden Berk et al.ll2004f T for the UV-optical variability. 

The comparison between the SF of galaxies, stars and 
QSOs in Fig. I6T71 shows that F var (QSO) > F var (star) > 
F var (galaxy). 

6.2. Cross match with GCVS 

The full samples of galaxies, stars and QSOs are 
matched with the General Ca talogue of Variable Stars 
(GCVS. iKholopov et al.l[l99l Vol. I-III and references 
therein) as an assessment to the definition of variabil- 
ity. The RA and DEC are both matched to within 1", 
respectively. 

Table[T2lshows the objects that exist in both the GCVS 
and our samples. As expected, all of the matched ob- 
jects are spectroscopically classified as stars by the SDSS. 
Three of the five GCVS stars are defined in this work to 
be variable candidates (> 3a detection), hence a com- 
pleteness of 60% when using 2-epoch observations to in- 
fer the stellar variability. 

The above indicates that the variability of stars and 
galaxies in this study could in principle be lim ited by the 
small number of epochs, also pointed out by llvezic et al.l 
(2003); and that the defining of variable stars by more 
than 3ct over the average in the variability may need to 
be adjusted. The fact that more than half of the vari- 
able candidates are real variable stars according to the 
GCVS shows that the calibration refinement is successful 
in selecting those objects. 

6.3. Cross match with SDSS Southern Stripe 

Our full samples of galaxies, stars and QSOs are also 
cross-matched wit h the variable catalog constructed by 
ISesar et all (|2007h based on the SDSS repeated imag- 
ing (the Stripe 82). The data cover more observation 
epochs (from 4 to 28, with an average of 9), useful for 
determining the completeness of our variable candidates 
based on 2 epochs. The matched objects are found to be 
eit her stars or QSOs in the SDSS spectral classification, 
as lSesar et al.l (|2007l )'s catalog is for un-resolved objects 
only. In Table [T3l the matched stars are listed. Five out 
of 22 sources matched to the SDSS variable catalog are 
defined to be variable stars in this work (> 3a detec- 
tion), implying a 23% completeness when using 2-epoch 
observations to infer the stellar variability. 

For QSOs, because more than 300 matches are found, 
only the number of found objects is summarized in Ta- 
ble [T4J along with our detection significance as variable 
candidates. About 10% of the objects are considered to 
be real variables in this work. If some of the narrow-line 
AGNs do vary in time, the completeness of the variable 
candidates from our current 2-epoch sample may be as 
low as 10%. 

7. SUMMARY 

We probe the spectroscopic variability in the galaxies 
and narrow-line AGNs in the optical wavelengths using 
multi-epoch observations in the SDSS. In order to detect 
the expected low amplitude of variability, we compare 
several approaches to refine the spectrophotometric cal- 
ibration in the SDSS. The final calibration, in terms of 



the mean and the standard deviation of the variability 
in our galaxy sample, is a factor of two better than the 
official SDSS reductions. 

Our sample of galaxy pairs spans in the restframe 
~ 110 days on average, with a maximum of ~ 700 days. 
The average change in the narrow-line AGN continuum 
flux, when converted to synthetic AB absolute mag- 
nitudes using the SDSS filters, is 0.07 ± 0.26 mag in 
the synthetic g band, where the uncertainty is prop- 
agated from the uncertainty in the spectral flux mea- 
surement. The fact that no obvious continuum variabil- 
ity is found is consis tent with the expectation from the 
AGN unified model (|Antonuccil 11993 : lUrrv fc Padovanil 
1995), in which some of the light-emitting (in our con- 
text, continuum-generating) structures in a narrow-line 
AGN, such as the accretion disk, are obscured by a dusty 
torus. Hence, any variability in such structures (cf. the 
UV-optical variability in QSOs in the time scale of 1 year) 
is not observed. This also provides an empirical evidence 
for modeling the continua of n arrow-line AGNs using 
stellar populati on models (e.g., Kauffm ann et al.l 120031 : 
lHao et al.ll2005h . Further, if we use this continuum vari- 
ability measurement to set an upper limit on the AGN 
activities, then the ratio between the flux of any varying 
spectral component, presumably related to AGN activi- 
ties, and that of the host galaxy is at most ~ 30%. 

By comparing the narrow emission-line ratios 
log 10 ([N II] /Ha) and log 10 ([O HI]/H/3) between the 
defined dim and bright phases in the continuum variable 
candidates, we find no evidence for their variability to 
be substantially larger than the la sample scatter in 
the line-ratio difference. This is consistent with the 
common view in which the narrow-line emissions in 
the AGNs are generated from low-density clouds. For 
example, an electron density n e ~ 1,000 cm~ 3 gives 
rise to a hydrogen recombination time ~ 100 years, 
which is typically much longer than the duration of the 
observation. 

We found no evidence of continuum variability in 
galaxies of eClass types A, B, C, D that is larger than 
the variability uncertainty. This agrees with the expec- 
tation that stellar light dominates their continuum, as in 
the narrow-line AGNs. 

We tabulate the upper limits of both the continuum 
and emission-line variability as a function of spectral 
type, using the classification schemes eClass and the BPT 
diagram. These values can serve as a sanity check for re- 
searchers who study variability for other sources. 

7.1. Variability in Classification 

There were few attempts in examining the galaxy 
classification from the poi nt of view of variabil- 
ity. A recent work by iBrunzendorf &: Meusingerl 
(2001) suggested variability is an efficient method to 
find narrow-line AGN. However, this promising re- 
sult was not subsequently con firmed by the authors 
(jMeusinger &: Brunzendora [2002). Our result on narrow 
emission-line ratios being non-variable seems to agree 
with the latter claim made by the authors. 

7.2. Next Steps 

So far we have focused on the average variability prop- 
erties of the galaxies and narrow-line AGNs. We are in- 



vestigating, on the object-to-object basis, any dramatic 
variation in spectral features or types. 
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Fig. 1. — The calibration spectrum (gray crosses) of the spectroscopic plate 390 and its 5th-ordcr polynomial fit (black line). The 
correction in flux density between two epochs is a function of the observed-frame wavelength. 
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Fig. 2. — An example star-forming galaxy without any calibration refinement (F var = 0.058). Top: spectra at the first (black) and the 
second (gray) epoch. Bottom: difference spectrum. 




Fig. 3. — The same galaxy as in Fig. [2] which shows zero continuum variability (-FVar = 0) after the SYS+BAND calibration refinement. 
Top: spectra at the first (black) and the second (gray) epoch. Bottom: difference spectrum. 
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Fig. 4. — -F V ar vs. S/N of galaxy spectra before (left) and after (right) the SYS+BAND calibration refinement. The number of objects 
in each spectroscopic plate is listed in Table [T] 
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Fig. 5. — The average variability of galaxy spectra over a given spectroscopic plate plotted versus the average spectral S/N. The spectra 
at S/N > 10 are selected in the subsequent definition of variable and non- variable candidates to minimize variation in F vaT among plates. 
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Fig. 6. — Restframe structure function of the QSOs, stars and galaxies. Only spectra with S/N > 10 are shown. Error bar is lcr m for 
each axis. Dotted line shows the best-fit structure function of an exponential functional form in Ar res t. 
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Fig. 7. — Restframe structure function. Same as Fig. [6] except that the fraction change in total flux |A///| is used instead of F v: 
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Fig. 8. — Restframe structure function of the QSOs, stars and galaxies. Different from Fig. [6] which shows the total variability, here the 
variability at the zero time lag is quadrature subtracted, respectively in each spectral type. Only spectra with S/N > 10 are shown. Error 
bar is la m for each axis. Dotted line shows the best-fit structure function of an exponential functional form in Ar res t. 
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Fig. 9. — The BPT diagram of the continuum variable candidates. The object type from the lowest to the highest [N II]/Ho value is: 
star-forming, composite and Seyfert 2. In each duo the filled square represents the average line ratios of the brighter objects, and the empty 
square for the dimmer ones. The solid and dashed lines respectively define the Seyfert 2 and star-forming galaxies as given by Kcwlcy ct al. 
(2006). The line ratios in all object types do not show substantial change between both epochs. The error bar represents the la sample 
scatter in each phase. 




Fig. 10. — Similar to Fig. [9] but in each duo the filled square represents the average line ratios of the first epoch, and the empty square 
for the latter epoch. The error bar represents the lc sample scatter in each epoch. 
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Fig. 11. — Similar to Fig. [9] but in each duo the filled square represents the average line ratios of the continuum variable candidates, and 
the empty square for the non-variable ones. Line ratios in both cases are from the bright phase. The error bar represents the 1<t sample 
scatter. 



TABLE 1 Spectroscopic plates used in the analysis of galaxies. 



plate a,b first epoch MJD second epoch MJD difference in MJD S 



266 


51602 


51630 


28 


a n t 

447 


356 


279 


51608 


51984 


376 


430 


361 


282 


51630 


51658 


28 


346 


278 


284 


51662 


51943 


281 


425 


336 


285 


51663 


51930 


267 


417 


317 


291 


51660 


51928 


268 


444 


344 


296 


51665 


51984 


319 


448 


364 


297 


51663 


51959 


296 


453 


369 


300 


51666 


51943 


277 


386 


288 


301 


51641 


51942 


301 


370 


282 


304 


51609 


51957 


348 


398 


317 


306 


51637 


51690 


53 


435 


352 


309 


51666 


51994 


328 


371 


287 


340 


51691 


51990 


299 


385 


302 


348 


51671 


51696 


25 


313 


240 


351 


51695 


51780 


85 


395 


311 


385 


51783 


51877 


94 


326 


259 


390 


51816 


51900 


84 


392 


313 


394 


51812 


51913 


101 


347 


267 


404 


51812 


51877 


65 


421 


335 


406 


51817 


52238 


421 


319 


239 


410 


51816 


51877 


61 


309 


243 


411 


51817 


51873 


56 


398 


318 


412 


51931 


52258 


327 


396 


306 


413 


51821 


51929 


108 


351 


274 


414 


51869 


51901 


32 


383 


299 


415 


51810 


51879 


69 


366 


277 


416 


51811 


51885 


74 


332 


259 


418 


51817 


51884 


67 


367 


285 


419 


51812 


51879 


67 


355 


293 


422 


51811 


51878 


67 


465 


375 


425 


51884 


51898 


14 


370 


290 


437 


51869 


51876 


7 


304 


235 


440 


51885 


51912 


27 


425 


331 


483 


51902 


51942 


40 


324 


247 


547 


51959 


52207 


248 


347 


288 


588 


52029 


52045 


16 


435 


344 


594 


52027 


52045 


18 


355 


285 


616 


52374 


52442 


68 


450 


349 


640 


52178 


52200 


22 


405 


323 


644 


52149 


52173 


24 


408 


339 






coi 70 

DZ1 1 O 


ol 


^too 


OOO 


712 


52179 


52199 


20 


216 


131 


790 


52433 


52441 


8 


402 


323 


803 


52264 


52318 


54 


430 


355 


804 


52266 


52286 


20 


422 


354 


810 


52326 


52672 


346 


502 


360 


814 


52370 


52443 


73 


363 


279 


820 


52433 


52438 


5 


307 


242 


960 


52425 


52466 


41 


427 


354 


972 


52428 


52435 


7 


417 


345 


978 


52431 


52441 


10 


328 


262 


1028 


52562 


52884 


322 


466 


371 


1034 


52525 


52813 


288 


471 


385 


1037 


52826 


52878 


52 


508 


417 


1291 


52735 


52738 


3 


421 


337 


1512 


53035 


53742 


707 


144 


99 


1782 


53299 


53383 


84 


457 


373 


2009 


53857 


53904 


47 


449 


354 


2351 


53772 


53786 


14 


332 


259 










23,330^ 


18,435 u 



Continued on next page 



There are 60 spectroscopic plates used, taken from the SDSS DR6. 

The plates 1664, 1905, 1907, 2075 (not listed in the table) are not used because the number of galaxy pairs is less than 
The total number of galaxy pairs in a plate. 
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d The number of galaxy pairs with successful calibration refinement using SYS+BAND. 
e The number of galaxy pairs summing over all plates, respectively. 
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TABLE 2 

Statistics of F var in galaxies using different calibration refinements. 



method 


mean 1 


0" m 


a 


No refinement 15 


0.020 


0.0002 


0.025 


SYS+BAND 


0.010 


0.0001 


0.009 


SYS+BAND (280 km S - 1 ) c 


0.012 


0.0001 


0.011 


SYS+OIII 


0.016 


0.0001 


0.015 


SYS+BAND (RED galaxies) d 


0.010 


0.0001 


0.009 


SYS+BAND (BLUE galaxics) d 


0.010 


0.0001 


0.009 



Note. — There are 11,766 common spectral pair in each of the calibration refinement approaches. Each spectrum is chosen to be refined 
successfully in each and every of the above procedures. In principal there can be more objects with successful refinement in each case. 



TABLE 3 

Upper percentile of the galaxy variable candidates' 1 . 



lower S/N 





5 


10 


15 


20 


upper S/N 


5 


10 


15 


20 


40 


numbcr b 


11 


102 


246 


206 


86 


percentile 


1.53 


1.89 


3.92 


6.26 


3.44 



Note. — In the galaxy sample using SYS+BAND refinement. 



TABLE 4 

Statistics of F var in the full galaxy sample. 



type 


number a,b 


mean 




<r a 


Seyfert 2 


92 


0.009 


0.0008 


0.008 


Starforming 


1333 


0.012 


0.0003 


0.011 


Composite 


314 


0.010 


0.0005 


0.009 


eClass A 


729 


0.010 


0.0004 


0.010 


eClass B 


919 


0.010 


0.0003 


0.010 


eClass C 


587 


0.011 


0.0004 


0.010 


eClass D 


483 


0.012 


0.0005 


0.010 


eClass E 


3 


0.024 


0.0100 


0.017 


eClass F 


5 


0.044 


0.0124 


0.028 


RED 


9603 


0.009 


0.0001 


0.009 


BLUE 


2705 


0.012 


0.0002 


0.011 



Note. — In the galaxy sample using SYS+BAND refinement. 



TABLE 5 

Statistics of F var in the galaxy variable candidates. 



type 


number a 


mean 


Urn 


<r a 


Seyfert 2 


5 


0.026 


0.0030 


0.007 


Starforming 


62 


0.037 


0.0023 


0.018 


Composite 


13 


0.031 


0.0043 


0.015 


eClass A 


34 


0.031 


0.0027 


0.016 


eClass B 


45 


0.031 


0.0020 


0.014 


eClass C 


37 


0.030 


0.0015 


0.009 


eClass D 


18 


0.034 


0.0030 


0.013 


RED 


388 


0.028 


0.0006 


0.012 


BLUE 


154 


0.034 


0.0012 


0.015 



Note. — In the galaxy sample using SYS+BAND refinement. 
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TABLE 6 

Difference in AB synthetic absolute magnitude between the dim and bright phases in g, r, i bands, for variable candidates. 



















< 


< 


< 




< 
< 


< 
< 


type 


(M g )* 


<M r ) 


(Mi) 


<*} b 


(AM,) 1 


{AM r ) 


(AMi) 


to 






b 


b 


b 


Scyfcrt 2 


-20.57 


-21.36 


-21.30 


0.14 


0.07 


0.05 


0.05 


0.26 


0.12 


0.20 


0.01 


0.03 


0.03 


Starforming 


-18.79 


-19.37 


-19.46 


0.09 


0.10 


0.05 


0.05 


0.09 


0.04 


0.08 


0.05 


0.03 


0.05 


Composite 


-19.79 


-20.53 


-20.53 


0.11 


0.09 


0.04 


0.03 


0.19 


0.08 


0.13 


0.05 


0.02 


0.03 


eClass A 


-20.70 


-21.55 


-21.25 


0.15 


0.08 


0.05 


0.06 


0.11 


0.05 


0.85 


0.04 


0.04 


0.05 


eClass B 


-20.12 


-20.95 


-20.74 


0.13 


0.09 


0.05 


0.04 


0.10 


0.04 


4.82 


0.03 


0.03 


0.03 


eClass C 


-19.95 


-20.73 


-20.81 


0.12 


0.09 


0.04 


0.04 


0.12 


0.05 


0.08 


0.03 


0.02 


0.03 


eClass D 


-18.18 


-18.75 


-18.91 


0.07 


0.09 


0.04 


0.06 


0.16 


0.07 


0.12 


0.05 


0.02 


0.06 


RED 


-20.14 


-20.99 


-20.91 


0.13 


0.08 


0.04 


0.05 


0.03 


0.01 


3.16 


0.04 


0.03 


0.04 


BLUE 


-19.27 


-19.87 


-19.89 


0.10 


0.09 


0.05 


0.05 


0.06 


0.03 


0.05 


0.05 


0.03 


0.05 



TABLE 7 

Best-fit characteristic variability time scale. 



type 


zero-point -F V ar a 


r a (year) b 


asymptotic F var (Ar res! - 


+ oo) reduced \' 2 


range in 


AT r[!at (day) 


QSO 


0.012 ±0.003 


0.4 ± 1.0 


0.022 ±0.002 


0.04 





- 300 


star 


0.010 ± 0.069 


11.2 ±6.3 


0.035 ± 0.035 


3.41 





- 710 


galaxy 


0.009 ±0.162 


8.5 ± 1.6 


0.025 ±0.081 


19.52 





- 500 



TABLE 8 

Best-fit characteristic variability time scale, after quadtrature subtraction of variability at zero time lag. 



type 


r s (ycar) a 


asymptotic F V ar(Ar rES i- 


+ oo) reduced x' 2 


range in 


Arrest (day) 


QSO 


0.4 ±0.0 


0.018 ±0.005 


13.10 





- 300 


star 


0.8 ±0.0 


0.012 ±0.003 


49.71 





- 710 


galaxy 


1.4 ±0.0 


0.011 ±0.018 


501.97 





- 500 



TABLE 9 

BPT line-ratio difference between dim and bright phases, for variable candidates. 



type 


(Alog 10 ([N II]/Ha)) a 


<Alog 10 ([O III]/H/3)> 


o-(Alog, n ([N II]/Ha)) b 


a(Alog 10 ([O III]/H/3)) 


Seyfert 2 


0.02 ±0.03 


0.06 ± 0.08 


0.03 


0.03 


Starforming 


0.09 ±0.01 


0.16 ±0.03 


0.17 


0.21 


Composite 


0.03 ±0.03 


0.11 ±0.09 


0.03 


0.08 


eClass C 


0.17 ±0.05 


0.47 ±0.12 


0.18 


0.43 


eClass D 


0.07 ±0.02 


0.15 ±0.06 


0.08 


0.17 



TABLE 10 
Statistics of F var in stars. 



method 


mean 




£7 


No refinement 


0.035 


0.0011 


0.053 


SYS±BAND 


0.012 


0.0003 


0.013 



Note. — See caption in Table [2] for the meaning of the statistics. There are 2,379 common spectral pairs in each of the refinement 
approaches. The statistics of F vaT are calculated at spectral S/N > 10. 



TABLE 11 
Statistics of F var in QSOs. 



method 


mean 


Cm 


a 


No refinement 


0.036 


0.0012 


0.034 


SYS±BAND 


0.015 


0.0005 


0.014 



Note. — See caption in Table [2] for the meaning of the statistics. There are 766 common spectral pairs in each of the refinement 
approaches. The statistics of F var are calculated at spectral S/N > 10. 



TABLE 12 

Matched objects between GCVS a and our samples. 



RA b (degree) 


DEC b (degree) 


spccClass c 


S/N 


Ar d 


difference in g 


difference in r 


difference in i 


-Fvar 


nSigma* 


218.750940 


-0.768439 


1 


20.4 


53 


0.13 ± 0.49 


-0.05 ± 0.30 


-0.09 ± 0.41 


0.017 


1 


118.540000 


42.817851 


1 


39.0 


7 


-0.02 ± 0.26 


0.00 ± 0.16 


0.00 ± 0.21 


0.005 





119.352760 


43.207637 


1 


43.0 


7 


0.05 ± 0.23 


0.01 ± 0.16 


-0.02 ± 0.22 


0.011 


3 


120.915900 


42.512477 


1 


46.5 


7 


-0.23 ± 0.20 


-0.03 ± 0.13 


0.08 ± 0.18 


0.051 


4 


133.434070 


57.811289 


1 


52.2 


40 


0.16 ± 0.18 


0.06 ± 0.12 


0.00 ± 0.18 


0.045 


4 



TABLE 13 

Matched objects between SDSS Stripe 82 variable catalog 3 - and stars in this work. 
c 

3 



Q ro S/N 3 At difference in g difference in r difference in i F VSLT nSigma 



356 


.032540 


1 


.027191 


1 


18.2 


10 


94 


0. 


03 


± 


1 


.64 





01 


± 


0. 


45 





.00 


± 





.49 





.003 





356 


.151320 





.303432 


1 


9.1 


12 


94 


-0. 


01 


± 


2 


.13 


-0. 


.02 


± 


0. 


68 


-0 


.04 


± 


0. 


73 





.007 





356 


.171110 


-1. 


220269 


1 


26.8 


10 


94 


0. 


08 


± 





.53 


0. 


.00 


± 


0. 


28 


-0 


.07 


± 





37 





.006 





354 


.635580 





.856453 


1 


7.9 


10 


94 





02 


± 


1 


.28 


-0 


.02 


± 





.94 


-0 


.00 


± 


1 


.68 





.039 


1 


12 


.348661 


-0. 


299891 


1 


52.6 


8 


101 





01 


± 





23 


0. 


.00 


± 


0. 


15 





.00 


± 





23 





.000 





13 


.369355 


-0 


722569 


1 


45.5 


8 


101 


-0 


13 


± 





.35 





.01 


± 


0. 


19 





.05 


± 





.27 





.025 


1 


31 


.675864 





.988745 


1 


17.0 


14 


65 


-0 


01 


± 





.60 





.00 


± 


0. 


35 





.01 


± 





.53 





.001 





36 


.891406 





.437566 


1 


11.7 


14 


421 





03 


± 





.79 





.03 


± 





.58 





.02 


± 





97 





.013 





35 


.787734 





.959517 


1 


25.6 


14 


421 


-0. 


06 


± 





.44 


-0 


.03 


± 





26 





.01 


± 





.37 





.015 


2 


36 


.347691 





.660534 


1 


20.4 


13 


421 


-0. 


04 


± 





.57 


-0 


.03 


± 





33 


-0 


.03 


± 


0. 


47 





.015 


2 


44 


.620576 


-0. 


609703 


1 


14.3 


14 


61 


0. 


18 


± 





.74 





01 


± 





45 


-0. 


.00 


± 


0. 


71 





.036 


4 


44 


.322878 





.785284 


1 


13.2 


14 


61 


-0. 


.14 


± 





.73 


-0 


.01 


± 





44 





.06 


± 





.68 





.031 


4 


46 


.057464 


-1. 


220639 


1 


29.6 


12 


56 


-0. 


.16 


± 





.32 


-0 


.03 


± 


0. 


20 





.07 


± 





.29 





.032 


4 


48 


.387960 





.715234 


1 


17.0 


14 


108 


-0. 


.07 


± 





.56 


-0 


.02 


± 


0. 


.35 





.03 


± 





.58 





.014 


1 


4!) 


.984166 


-0 


044137 


1 


8.6 


14 


108 


-0. 


12 


± 


1 


.16 


-0 


.01 


± 


0. 


.79 





.05 


± 


1 


39 





.030 





51 


.883714 





.064299 


1 


54.1 


12 


32 


-0. 


40 


± 





.19 


-0 


.05 


± 


0. 


12 





.15 


± 





.17 





.085 


4 


52 


.295747 





.603985 


1 


20.9 


12 


32 


-0. 


.08 


± 





.62 


-0 


.01 


± 


0. 


34 





.04 


± 





.49 





.015 


2 


52 


.476904 





.389118 


1 


45.5 


14 


32 





02 


± 





.28 





.00 


± 





.17 





.00 


± 





24 





.006 





51 


.602434 





.674580 


1 


45.2 


14 


20 





01 


± 





.24 





.01 


± 


0. 


.16 


-0 


.04 


± 





21 





.007 


4 


43 


.924804 


-0 


.542138 


1 


2.8 


14 


707 


-0. 


.10 


± 


4 


.41 





.15 


± 


2 


26 





.17 


± 


3 


21 





.017 





44 


.201206 





.013831 


1 


19.2 


14 


707 


-0. 


.15 


± 


1 


.10 


-0 


.01 


± 





30 





.03 


± 





30 





.011 





44 


.918114 





.758511 


1 


6.3 


12 


707 


-0 


14 


± 


2 


.89 


-0 


.17 


± 





98 


-0 


.28 


± 


1 


13 





.100 


2 



NOTE. — See also captions in Table [T21 



TABLE 14 

Matched objects between SDSS Stripe 82 variable catalog 3 and QSOs in this work 



nSigma b 


number of objects 


percentage 





231 


67.3 


1 


45 


13.1 


2 


33 


9.6 


3 


10 


2.9 


4 


24 


7.0 



